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Abstract 

This letter presents results of new high resolution = 1 Cold + Hot Dark Matter 
(CHDM) and Cold Dark Matter (CDM) simulations. Properties of groups in these sim- 
ulations reflect the lower small-scale velocities and the greater tendency to form distinct 
fllaments on both small and large scales in CHDM as compared to CDM. The fraction of 
galaxies in groups and the median group rms velocity are found to be powerful discrimina- 
tors between models. We combine these two features into a very robust statistic, median 
group rms velocity t^gr(/gr) as a function of the fraction fgj- of galaxies in groups. Using 
this statistic, we compare "observed" simulations to CfA data in redshift space in a careful 
and consistent way. We flnd that CHDM remains a promising model, with for example 
Ugr(0.45) ~ 125 lb 25 km in agreement with the CfA data, while CDM with bias b=1.0 
(COBE-compatible) or b=1.5, both giving Ug].(0.45) ~ 400 it 25 km s~"^, can be virtually 
ruled out. Using median M/L, the observed value of O is 0.10 (CHDM) to 0.38 (CDM). 

Subject Headings: cosmology: theory - dark matter - galaxies: clustering 
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1. Introduction 



A simple inflationary = 1 universe with a primordial Zeldovich spectrum of Gaus- 
sian, adiabatic fluctuations and a mixture of Cold + Hot Dark Matter (CHDM) has proven 
to be a serious candidate for the new, standard cosmological model (Klypin et al. 1993, 
hereafter KHPR), now that standard CDM seems unable to account simultaneously for 
small-scale galaxy pairwise velocities, large-scale motions and structure, and observed CMB 
anisotropies. Following KHPR, we consider a model with Ocoid = 0.6, Oj/ = 0.3, Obaryons = 
0.1. The hot component is assumed to be a light neutrino, and specifying the model's 
single additional free parameter, the hot dark matter density fraction Oj^, flxes the neu- 
trino mass. Interestingly, in the MSW scheme, the simplest and most attractive solution 
of the solar neutrino problem, the muon neutrino mass must be about 3 milli-eV; the 7 
eV mass implied by Oj^ = 0.3 is then in the range predicted for the tau neutrino by the 
simple (albeit speculative) "seesaw" model relating the masses of the neutrinos to those of 
the quarks (see e.g. Ellis eA al. 1992). 

When the CHDM initial fluctuation spectrum is normalized to the COBE Qps-norm = 
17/iK (Smoot et al. 1992), corresponding to linear biasing parameter h = 1.5 (erg = = 
0.67), the CHDM model is found to compare favorably to observations on a number of mea- 
sures, much better overall than CDM. These include small scale galaxy velocity dispersion 
and large scale streaming velocities (KHPR), cluster-cluster correlation function (Holtz- 
man and Primack 1993, Klypin and Rhee 1993), power spectrum from CfA2 (Vogeley eA 
al. 1992) and IRAS (Fisher et al. 1993, Feldman et al. 1993, Klypin et al. 1993b). 

KHPR describe results of the analysis of particle-mesh (PM) simulations on a 256^ 
mesh for 14, 50, and 200 Mpc boxes (Hubble parameter ifo = 50 km Mpc ~^ is assumed 
throughout this paper). Klypin et al. (1993b) describes a more complete analysis of new, 
higher resolution CHDM and CDM simulations. Here, we focus on the most striking 
results from our study of properties of galaxy groups in these high resolution simulations, 
and compare them to identically selected groups in CDM simulations and in those formed 
from the CfAl Survey (Davis et al. 1982). 

2. Simulations, Galaxy Identification 

The CHDM simulations were evolved from the initial power spectrum appropriate 
for the CHDM model at 2; = 15 to the present. Calculations were done with a PM code 
run on a 512^ mesh in a 100 Mpc box (i.e. cell size = 195 kpc) with 256^ cold particles 
(Mparticie = 2.9 X 10^ M©), and 2 x 256^ hot particles (Mparticie = 6.2 x 10^ M©). For 
comparison, we also ran simulations of CDM (same mesh, Mparticie = 4.1xlO^M0 ) at the 
same bias factor h = 1.5 (hereafter CDM1.5), and at the COBE-compatible CDM bias 
of 1.0 (hereafter CDMl). CDM1.5 was begun at z = 18, CDMl at z = 27.5. Due to a 
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statistical fluke of probability about 1/10, our flrst 100 Mpc CHDM simulation (hereafter, 
CHDMi) had more power on the 100 Mpc wavelength than would be typical for a 100 
Mpc box. The CDM simulations were run from the same choice of input waves. We also 
report here results from another simulation, CHDM2, which had a more typical power 
spectrum for this box size. After running these simulations, we found that there were two 
mistakes in our initial conditions (see KHPR, Note Added in Proof): the fltting formula 
for the cold spectrum was too small and the velocities were too large, both by about 20% 
on small scales. However, these effects are in phase and largely cancel. We have run 
another CHDM2 (rev) simulation with both errors corrected, and found that the power 
and velocity differences between this and CHDM2 (old) on small scales declined to 5% by 
z = 1 and remained at this level. All results discussed here are from CHDM2 (rev), but 
the argeement with CHDM2 (old) is well within the la error bars; thus the CHDMi and 
CDM simulations should also be reliable. 

Dark halos ("galaxies") in our simulations were identifled as local density maxima 
on the 512^ mesh, and assigned the mass inside this cell. Only dark halos with mass 
M > 4.1 X 10"^" Mq, corresponding to Sp/p > 80, were used in the analysis below. 

3. Sky Catalogs, Group Identification 

Group properties reflect features of the CHDM and CDM models on the ~ 0.1 — 10 
Mpc scales where they differ most signiflcantly. In order to make reliable comparisons with 
observed groups, it is essential to "observe" simulations and identify groups in redshift 
space in exactly the same way as is done for the data. Since our simulation resolution 
of 195 kpc per cell is more coarse than observations, we began by merging CfAl grouped 
or paired galaxies from Nolthenius (1993; N93) which were closer than a sky projected 
separation of 235 kpc. This merged catalog had 2112 galaxies, compared to the original 
2406. 195 kpc is only ~ 10% of the median group harmonic radius (r/j)nied, and our limited 
resolution should not seriously affect our results. (See NW for deflnitions of {rh)ined and 
other quantities used here.) 

We assigned luminosities L to halo masses in two different ways. The flrst and in 
our opinion more reliable prescription was to randomly generate Schechter-distributed 
magnitudes, using merged CfAl parameters a = —1.17 and M* = —21.25, and assign these 
to halos with luminosity monotonically increasing with halo mass (Schechter method). We 
also tried using the relation between dark matter density and baryonic density found in 
hydrodynamic simulations by Cen and Ostriker (1992, 1993), assuming constant Mbaryon/-^ 
(CO method). (We used the masses in the 3^ cells centered on each halo for p = p/p in 
CO's relation log(pbaryon) = A + 2.3 log(ptot) -0.20 log^(ptot), and took A = 5.8 for CDM 
and A = 6.23 for CHDM simulations to match the merged CfAl galaxy density.) The CO 
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method produces a luminosity function which is too steep at low luminosities (a ~ —1.8, 
essentially the Press- Schechter result) and has too many bright galaxies. This weighted 
the sample in favor of distant, bright (and to a lesser extent, nearby and faint) galaxies 
much more heavily than the CfA data, and moderately raised median group sizes, rms 
velocities, and fractions grouped. While we regard the CO L{M) as unrealistic, it turned 
out to give essentially the same results for our favored statistic t^gr(/gr) as our preferred 
Schechter method. 

We partially compensate for cosmic variance by selecting six viewing locations which 
approximate our place in the local universe, using the same viewing location criteria 
for both CHDM and CDM simulations: (a) after applying a magnitude limit (using the 
Schechter method), halo density within 15 Mpc in redshift space was within a factor 1.5 
of CfAl density, and (b) a roughly Virgo-mass small cluster was ~ 25 Mpc distant. (By 
using the same random numbers for the CHDMi and CDM simulations, we insured that 
the relative halo densities and locations of large concentrations were virtually the same, 
and we could then place CDM viewers on the nearest halo to the CHDMi viewers.) As- 
signed magnitudes were then cut at m = 14.5, and the resulting catalogs randomly culled 
an additional ~ 20% down to CfAl density. Finally, in order to include edge effects, a 
20 degree wide "zone of avoidance" was added. This produced sky catalogs of 10.38 sr, 
compared to the 2.66 sr of the CfAl sample. 

We then used the adaptive friends-of- friends link algorithm of N93 to identify groups. 
Each galaxy is surrounded by a cylinder with its axis along the line of sight, and three or 
more touching or intersecting cylinders are considered to define "groups." The radius of 
the cylinder is a constant fraction Dn ("sky link" parameter) of the mean galaxy spacing on 
the sky, and the length is VL{d) ("redshift link," parameterized by its value V5 = ¥^(5000) 
at the typical CfAl distance d = 5000 km s~"^). Both sky link and redshift link increase 
with sample incompleteness (and therefore distance) according to the fitted luminosity 
functions for each catalog as discussed in N93. A fiducial Dn = 0.36 (1.2 Mpc at distance 
1000 km s~^) produces groups of minimum overdensity about 20 and is a good compromise 
between selecting only dense cores or percolating connections between groups on the sky; 
it is quite similar to that used in earlier studies (Geller and Huchra 1983; Nolthenius and 
White 1987 (NW); RameUa, Geller, & Huchra 1989; N93; and Moore et al. 1993 (MFW); 
note that NW, MFW, and N93 express Dn as a fraction of the 3-D particle spacing, rather 
than the more appropriate projected spacing used here, so MFW's fiducial Dq = 0.3 
corresponds to our Dn = 0.3 X 7r/2 = 0.47). 
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4. Results 



The primary results of this study foUow from two findings: CHDM simulations are 
both much more filamentary and also have lower small-scale velocities than the CDM 
simulations. This means that groups form and percolate much more readily in CHDM, both 
in real and redshift space. Nevertheless, we find that the full distributions of multiplicity, 
mass, velocity dispersion, and size near the fiducial links are poor discriminators between 
models; all our simulations are similar, and, within the errors, in reasonable agreement 
with the CfAl data. 

NW first demonstrated the value of comparing the median/percolation properties of 
groups versus link parameters, finding that CfAl groups form and percolate more readily 
than in low density b=1.0 CDM. Here we extend this idea and present a sensitive and robust 
statistic for discriminating between model simulations. Note that for most purposes our 
conclusions do not require that "groups" actually be dynamically isolated objects, only 
that they be consistently identified. We may therefore make valid comparisons througout 
the link range, even though conventionally defined groups are only approximated near the 
fiducial links. In this sense "groups" are only a label we apply to a set of objects which 
provide useful discriminating statistics. 

Figure 1 compares the fraction of galaxies fgj- found in groups vs. redshift link V5 
at our fiducial sky link = 0.36. /gr rises steeply at small V5 as valid group members 
are added to groups. Within this regime, median group rms velocities Vgj- are artificially 
clipped to that allowed by V5; we refer to this as the "clipped regime". At higher V5, 
only interlopers from the lower density surroundings are picked up, and the curves flatten 
("interloper regime"). This transition is rather distinct for the CHDM and CfA data. This 
is consistent with our visual inspection of the halo distribution: we found that CHDM 
halos are mostly found in fllaments and knots, much more so than in the CDM boxes, 
which show a more gradual transition to the interloper regime. There is remarkably little 
variation between the six different viewpoints, as shown on the simulation points on all 
the flgures by the la error bars. For the CfA points, the errors were estimated from the 
CHDMi simulation; we took the rms spread among its six viewpoints after restricting 
each viewpoint's sky to the same 2.66 sr divided into two regions as for CfAl, one region 
including the nearby Virgo-sized cluster. Such error estimates include both statistical 
scatter and cosmic variance, but the difference between the two CHDM simulations shows 
that these error bars do not include all the cosmic variance. It is clear that both CDMl 
and CDMl. 5 seriously conflict with the data. 

Next we consider how the median group rms velocity Vgj- depends on the redshift link. 
(We deflne the rms velocity of each group as Vgj- = [(n — < ^ >)'^]^^'^ = 

"group velocity dispersion" in Geller and Huchra 1983; here Vi = czi.) Figure 2 shows 
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function of V5, here done at the larger = 0.47 to facihtate comparison with 
MFW. Similar to Figure 1, all curves show a steeper rise in the clipped regime followed 
by a shaller rise in the lower density interloper regime. A similar curve at our fiducial 
Dn = 0.36 shows CHDM transitions at V5 ~ 350 km s'^ and CDM at V5 ~ 600 km s'^ 
The curves show that within the clipped regime group rms velocities are limited by the V5 
link rather than the underlying properties of the simulation. The transition to a flatter 
slope is slower in the CDM curves of Figures 1 and 2, indicating a more substantial tail of 
high velocity neighbors. CDM Ug^'s are signiflcantly and consistently higher than the data 
even at low V5. 

For our conclusions, it is important to be sure we have not overestimated CDM small 
scale velocities. MFW used the CDM simulations of Frenk et al. (1990) to study galaxy 
groups. Their P^M code gave better spatial resolution than our PM code, and their galaxy 
identiflcation scheme ("high peaks"), galaxy luminosity function, and link algorithm (NW) 
also differ from ours. Nevertheless, our Ugj. for CDM1.5 at the same Dn compares well with 
their CDM b=1.6 results, as shown in Figure 2. (No error bars are shown on the MFW 
points, since MFW did not estimate them.) 

Figures 1 and 2 suggest that the optimal V5 for producing the largest relatively un- 
contaminated groups is near V5 ~ 350 km s~^ for CHDM groups and 600 km s~^ for CDM 
groups. To conflrm this, we made groups from these same sky catalogs using full 3D infor- 
mation with a 3D link parameter Dn = 0.36. Reproducing the 3D group's Ug^'s required 
a 1/5 ~ 350 km s'^ for both CHDM simulations, and a ^5 = 450 km s'^ (CDM1.5) to 
550 km (CDMl) for CDM. Matching fractions grouped required slightly higher Vs's. 
We cover this range by adopting flducial V5 of 350 for CHDM and 600 km s~^ for CDM. 
Note that V5 = 350 km also gives the optimum CfAl groups (NW, N93). 

A statistic that is both discriminatory and robust is obtained by plotting the median 
group rms velocity as a function of the fraction of galaxies in groups. Since raising V5 
typically raises both /gj. and Ugr, plotting t^gr(/gr)|D„ at flxed Dn (i.e. eliminating V5) 
produces a strong discriminator between models, shown in Figure 3. From left to right. 
Figure 3 points correspond to V5 from 40 to 2000 km . Note that through its implicit 
dependence on V5, this plot is sensitive to differences in the small-scale velocity dispersion 
between models. This statistic turns out to be quite robust. Using the very different Cen- 
Ostriker luminosity assignment produces virtually identical curves, as does a change from 
Dn = 0.36 to 0.47 (though here all curves are pushed to slighly higher /gr)- It does not 
appear possible to reconcile CDM groups with the CfA data for any reasonable selection 
parameters or mass-to-light assignment. Over most of the range, CHDM agrees well with 
the CfA data. Note also that, to the extent the CfA region is typical, the CHDM2 initial 
conditions are a more appropriate comparison. On all flgures, CHDMi (and associated 
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CDM curves as a block) should then be moved towards the CHDM2 curves. This worsens 
CDM's and improves CHDM's agreement with the data. 

We also considered a related statistic t^gr(/gr)| V5 , this time varying the sky link pa- 
rameter Dn, keeping V5 = 600 km fixed (see Klypin et al. 1993b). Again, CDM Ug^'s 
are 50 — 100 km too high and clearly inconsistent with the data. CHDM curves agree 
much better, but are moderately too high at small /gj. (corresponding to small Dn)- 

Finally, we determined the observed value of O from the simulations. Using the usual 
MIL method (e.g. NW), we found O's of 0.10 (both CHDM), 0.17 (merged CfAl), 0.30 
(CDM1.5) and 0.38 (CDMl). CDM O's are higher largely due to their higher fiducial 
V5 = 600 km . These demonstrate that it is indeed possible to measure low O from 
groups in an O = 1 universe. 

5. Caveats and Conclusions 

While CHDM appears to do much better than CDM in these comparisons to obser- 
vations, there remain discrepancies that appear significant. In particular, at small Dn 
we are looking at dense group cores, which appear to be significantly cooler in the CfA 
data than any of the simulations. At the smallest Dn = 0.15, CHDM Ug^'s are almost 
50 km higher than the data. Also, CfA Ug^'s show a significantly steeper rise with 
/gr than either CHDM or CDM. A related point is that simulation groups show larger 
median group harmonic radius (r/j)nied than those observed, by ~ 20% for CHDM, ~ 30% 
for CDM, especially at small V5 and Dn- These facts suggest that real group cores are 
both colder and denser than in CHDM and especially CDM simulations. The additional 
10% discrepancy seen in CDM is likely to be genuine, a result of higher pairwise veloci- 
ties infiating the groups. For CHDM, the discrepancy may be due in part to our limited 
resolution at these scales {Dn = 0.15 ~ 0.75Mpc ~ 3 — 4 cells for nearby groups) or the 
fact that dissipation is not modelled, or may indicate a real difference. At least some of 
the (r/j)nied discrepancy is likely due to the problem of the "overmerging" of dark matter 
halos in dissipationless simulations (e.g. Gelb and Bertschinger 1993), which is undoubt- 
edly present in our simulations. Overmerging will tend to remove close neighbors and thus 
infiate average pairwise spacings. In CDM models this appears significant for halo masses 
M > 10^^ Mq (Evrard, Summers, & Davis 1993), which comprise 40% of the halos in 
our CDM models, most of these in groups. However, overmerging actually lowers pairwise 
velocities, so any correction would most likely worsen CDM's already desperate plight with 
respect to group rms velocities. For CHDM, the hot component delays collapse epochs and 
thus helps reduce overmerging, and while no studies comparable to Evrard, Summers, & 
Davis suggest at what mass it becomes significant, only 21 — 28% of our CHDM halos 
are above 10"^^ Mq. Thus, correcting for overmerging is likely to widen the differences 



7 



between CDM and CHDM, and between CDM and the data, thereby strengthening our 
basic conclusions. 

To summarize our conclusions: Median group rms velocity Vgj- versus fraction of galax- 
ies in groups fgj- is shown to be a powerful and robust statistic for discriminating between 
cosmological models. Both unbiased and low bias = 1 CDM groups show unaccept- 
ably high median rms velocities and too few galaxies in groups, due to the combined 
effects of high pairwise velocities and less tendency for halos to be arrayed in filaments 
and knots. CHDM's highly filamentary nature percolates more easily and makes groups 
slightly easier than observed in the CfAl data. CHDM remains a promising model, with 
Ugr(0.45) ~ 125 lb 25 km at = 0.36 in agreement with the CfAl data. At our fiducial 
Dn and V5, /gr = 0.42 for merged CfAl and there are 170 groups. A bootstrap estimate of 
the purely statistical error in Ug].(0.45) for CfAl is 9%; the rms among the six viewpoints 
is 11%, giving Ugr(0.45) = 138 ± 15 km s'^ CDM with bias b=1.0 (COBE-compatible) 
and b=1.5, both with Ug].(0.45) ~ 400 it 25 km at = 0.36, can be virtually ruled 
out. 

The Vgj-(fgj-) statistic should be measured for the larger redshift data sets that will 
soon become available (e.g. CfA2, SSRS2, Perseus-Pisces), and also on high-resolution 
simulations of all the potentially viable cosmological models, including CHDM with varying 
fractions of Oj^ 0.2 — 0.25, and Obaryon ~ 0.05. It will be interesting to see which models 
can pass this difficult test! 

Acknowledgments. We thank J. Ostriker for suggesting that we use CO for assigning 
luminosities, and the Aspen Center for Physics for its hospitality in August 1993. Our 
simulations were done on the Convex-3880 supercomputer at the NSCA. 
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U vs. V5 
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Figure 1. The fraction of galaxies grouped rises sharply at low redshift link V5 as valid group 
members are added, then levels off as only interlopers are included. CDM consistently 
groups a lower fraction than the data, and shows a less distinct transition to the interloper 
regime at higher V5. CHDM groups slightly too high a fraction. 
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Vgr vs. V5 
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Figure 2. is here plotted vs. V5 at MFW's fiducial with MFW's CDM results also 
plotted. Vgj- rises more steeply at low V5, then flattens somewhat for the same reasons as in 
Figure 1. At low V5, rms velocities are artiflcially clipped for all simulations. Our flducial 
Vs's (at our flducial = 0.36) are marked with arrows. Beyond the clipped regime, CDM 
rms velocities are consistently too high. 
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Vgr vs. /gr 




Figure 3. Median group rms velocity t^gr(/gr)|D„ of galaxies in groups at fixed sky link 
Dn proves to be a powerful and robust statistic. It is plotted here for our fiducial sky 
link Dn = 0.36, and Schechter luminosity assignments. CDM simulations for both bias 
choices produce Vgj- far higher than CfAl groups, while both CHDM simulations are in 
good agreement. Relative positions of curves remain virtually unchanged under changes 
in sky link or luminosity assignment method. 
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